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ABSTRACT 


Water  jets  are  produced  by  vertically  accelerating  a 
rotating  cone  partially  filled  with  water.  It  is  shown  that  the 
acceleration  of  the  parabolic  meniscus  results  in  amotion  sim¬ 
ilar  to  that  observed  in  a  shaped  explosive  charge.  (Monroe 
jet).  Acceleration  of  the  cone  is  effected  by  means  of  an  induc¬ 
tive  electromagnetic  accelerating  device  (conical  pinch)  whose 
theory  i3  developed  in  terms  of  the  WKB  approximation.  A 
second  order  inviscid  theory  is  presented  of  the  motion  of  the 
fluid  in  the  cone  in  terms  of  the  Penney-Price  linearization 
procedure  and  it  is  shown  that  good  agreement  for  the  jet  head 
velocity  can  be  achieved  for  low  velocities. higher  veloc¬ 
ities,  experimental  results  appear  to  lag behind  thfeory,  probably 
owing  to  the  dispersal  of  the  jet  head  through  viscous  drag  with 
the  surrounding  atmosphere.  The  shape  of  the  jet  at  early  times 
is  well  represented  by  first  order  theory. 
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HIGH  SPEED  WATER  JETS  FROM 
VERTICALLY  ACCELERATED  ROTATING  CONES 


1.0  INTRODUCTION 

High  speed  water  jets  have  for  some  time  been  used  for  drilling,  cutting  and 
shaping  various  materials.  The  advantages  of  this  technique  lie  in  the  small  width  of 
attainable  cuts  and  the  avoidance  of  the  necessity  to  sharpen  the  work  tool.  Since  most 
methods  of  water  jet  production  comprise  of  high-pressure  pump  followed  by  pipes,  noz¬ 
zles  and  possibly  hydraulic  intensifiers,  the  maximum  stagnation  pressures  attainable  in 
the  jets  depends  on  the  strength  of  the  walls  used  to  enclose  the  hydraulic  circuit.  To 
overcome  this  limitation,  it  has  been  suggested  that  some  form  of  hydrodynamic  jet  inten¬ 
sification  involving  converging  flow  may  be  used.  One  such  device  is  the  Monroe  jet  or 
shaped  charge.  In  this  case  an  explosive  charge  containing  a  conical  hole  lined  with  a 
layer  of  metal  is  exploded  at  the  end  nearest  the  apex  of  the  cone.  The  resulting  detona¬ 
tion  wave  compresses  the  metal  lining,  causing  it  to  flow  towards  the  centre  line  and  along 
it  to  form  a  narrow  jet  of  considerable  speed  and  penetrating  power.  The  theory  of  this 
mechanism  has  been  fully  explored  and  tested  experimentally ' '4l .  As  the  jet  is  free¬ 
standing  it  does  not  require  confining  walls  and  hence  the  problem  of  wall  strength  does 
not  arise.  Obviously  the  repeatability  of  such  a  device  is  severely  limited  by  the  fact  that 
the  source  of  the  jet  is  totally  destroyed  in  the  process.  For  this  reason  the  suggestion 
has  been  made  that  the  explosive  charge  be  replaced  by  a  magnetic  compression  device 
and  the  metal  lining  by  some  sort  of  liquid  film.  That  magnetic  pressures  of  great  mag¬ 
nitude  can  be  generated  has  recently  been  shown  by  Bless51.  It  remains  to  couple  the 
magnetic  force  in  some  way  to  the  liquid  film.  Preliminary  experiments  in  this  labora¬ 
tory,  using  a  layer  of  mercury  indicate  that  the  liquid  becomes  unstable  under  an  external 
magnetic  field  and  no  jet  is  formed.  Further  experiments  using  non-conducted  fluids 
contained  in  a  deformable  hollow  metal  cone  were  equally  unsuccessful,  as  the  cone  be¬ 
came  subject  to  elastic  instability,  except  in  the  case  of  small  strain,  as  will  be  further 
discussed  below. 

As  the  theory  of  the  Monroe  jet  uses  nothing  but  ideal  incompressible  and  in- 
viscid  fluid  characteristics,  it  may  be  argued  that  a  cone  lined  with  a  liquid  layer  and 
accelerated  along  its  axis  is  in  fact  equivalent  to  a  cone  that  collapses  in  a  manner  in 
which  all  points  on  the  cone  wall  approach  the  centre  line  with  equal  velocity,  the  only 
difference  being  that  in  the  former  case  an  additional  shear  velocity  along  the  cone  wall 
is  imposed,  which,  since  the  liquid  is  assumed  inviscid,  does  not  alter  the  mechanics  of 
the  motion.  This  is  the  solution  adopted  in  this  work,  whereby  the  liquid  layer  is  pro¬ 
duced  by  rotating  the  cone  about  its  axis.  It  will  be  shown  in  the  following  that,  although 
the  liquid  layer  produced  in  this  way  is  by  no  means  of  uniform  thickness,  the  velocity 
amplification  of  the  resulting  jet  relative  to  the  cone  is  considerable  and  in  the  case  of  the 
meniscal  parabola  being  tangent  to  the  cone,  corresponds  to  the  findings  of  the  Monroe 
jet  theory. 

2.0  ELECTROMAGNETIC  ACCELERATION  OF  A  CONE 

The  theory  of  the  "electromagnetic  cannon"  or  accelerator  has  been  extensively 
considered6  /),  among  others  by  Mostov  et.al.71  who  present  probably  the  most  complete 
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and  elegant  treatment.  We  shall  here  reproduce  some  of  the  basic  equations  but  then 
resort  to  an  analytical  approximation  involving  the  WKB  scheme,  rather  than  consider 
asymptotic  and  power-series  solutions  as  in  Reference  7.  The  analytical  foundations  of 
the  WKB  method  have  been  well  documented81  and  we  shall  show  by  means  of  comparison 
with  analog  solutions  that  it  yields  results  of  adequate  accuracy  in  the  present  case  also. 
We  depart  from  the  treatment  in  Reference  7  not  only  in  the  method  but  also  in  the  fact 
that  for  practical  reasons  we  have  to  take  account  of  the  fact  that  the  inductance  of  the 
circuit  is  not  a  linear  function  of  the  displacement  of  the  load. 


We  consider  at  first  a  circuit  made  up  of  a  condenser  of  capacity  c  charged  to 
a  voltage  V0  followed  by  a  circuit  resistance  R  and  in  inductance  L.  The  metallic  cone 
is  thought  to  be  imbedded  in  the  coil  whose  inductance  is  a  linear  function  of  the  axial 
displacement  x  of  the  cone,  i.e. 


L  =  L0  +  x.dL/dx  =  L0  +  xL, 


(1) 


The  quantity  of  greatest  interest  will  be  the  transfer  efficiency  e,  which  is  the  ratio  of  the 
kinetic  energy  imparted  to  the  cone  of  mass  m  to  the  energy  initially  stored  in  the  con¬ 
denser.  It  can  be  shown  that  the  magnetic  force  F8  driving  the  cone  is  given  by: 


Fg  -  L,i2/2 


where  i  is  the  instantaneous  current  in  the  circuit.  Hence,  we  have  the  circuit  equation: 


di 

dt2 


(Li) 


4 


i/c  =  0 


(2) 


and  the  force  balance  equation: 


m.jf  -  L.lVS 


(3) 


The  relations  (l)-(3)  form  a  system  of  equations  which  completely  describe  the  motion 
provided  a  suitable  set  of  initial  conditions  is  imposed.  To  find  the  transfer  efficiency 
we  shall  however  employ  the  energy  relation  directly.  If  the  velocity  of  the  cone  is  de¬ 
noted  by  v,  we  have  from  energy  balance: 


v’c/2  =  R  / i2dt  +  mv2/2 
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indicating  balance  between  condenser  energy,  kinetic  energy  of  the  cone  and  ohmic  heat 
dissipation.  Integrating  (3),  we  obtain: 


m.v  =  1/2L ,  J  i2dt 


Eliminating  the  integral  between  these  two  relations,  we  obtain  a  quadratic  relation  in  v 
Solving  in  terms  of  the  transfer  efficiency: 


e 


mv2/cVj 


2mR2 

EL,2 


[l  -  (1+  EL*/2mR2)H2  ]J 


(4) 


where  E  -  V„ c/2  is  the  condenser  energy.  In  this  treatment  it  is  assumed  that  the  cir¬ 
cuit  resistance  is  independent  of  the  position  of  the  cone,  an  assumption  which  is  usually 
borne  out  in  practice,  where  the  external  losses  in  the  conductors  and  the  condenser  die¬ 
lectric  predominate  over  the  losses  in  the  cone.  An  even  simpler  expression  for  the 
transfer  efficiency  can  be  obtained  by  introducing  a  velocity  characteristic  of  the  electro 
mechanical  circuit,  the  "circuit  speed": 


C  =  R/2L, 


(5) 


Though  this  expression  is  only  a  symbolic  velocity,  it  can  be  shown  that  by  means 
the  equation  governing  the  motion  of  the  driven  gas  in  an  eiectromagnetic  shock  tube  ca 
be  rendered  in  a  form  exactly  analogous  to  that  of  the  gas-driven  shock  tube  (Taub  s 
equation),  where  C  then  corresponds  to  the  sound  velocity  of  the  undisturbed  gas.  Corre 
spondingly,  one  can  define  a  '  circuit  Mach  number  : 


M  =  v/C 


in  terms  of  which  the  transfer  efficiency  becomes: 

e  =  M/(M  +  8) 


Within  the  assumptions  made,  Equations  (4)  and  (6)  are  exact;  however  for  the 
sake  of  the  generalization  made  in  the  following  it  will  be  useful  to  attempt  an  approximate 
treatment  tased  on  the  WKB  scheme.  This  will  enable  us  not  only  to  calculate  the  trans¬ 
fer  efficiency  for  the  more  general  case  of  the  Inductance  not  being  a  linear  function  of 
the  cone  displacement,  but  also  to  study  the  motion  and  the  development  of  the  electrical 
quantities  with  the  passage  of  time.  It  is  convenient  to  first  transform  the  relevant  equa¬ 
tions  into  non-dimensional  units,  thus: 


T  =  t/(cL0)'12 ,  Z  =  xL,/L0,  K  -  R.  (c/L0)’12  *  Y  iL,(c/2m) 


112 


(7) 
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Eliminating  L  between  (1)  and  (2),  one  can  now  write  (2)  and  (3)  in  the  form: 


iP  V(1  *  zi  *  K  g  +  Y  =  0 


(8) 


O.  _  2 

,2  -  I 


dT' 


(9) 


The  following  initial  and  boundary  conditions  apply: 


Y  =  0 


Y'  =  P 


Z  =  0 


r  when 


Z'  =  0 


T  =  0 


(10) 


Y  -*  0  as  T  — ’  *> 


where  primes  denote  differentiation  with  respect  to  T.  From  the  electrical  conditions  at 
T  =  0  it  follows  that: 


P  =  V0cL,/(2mL0)''2 


(11) 


To  eliminate  the  first-order  term  in  (8)  we  set: 


Y  =  (1  +  Z)*,.g(D.exp.(-  ^ 


where  g(T)  is  an  unknown  function  of  T.  One  obtains  from  (8): 


i  K71 

g"  +  g(i  +  Z)-2  (1  +  z-  Y~  -  —)  =  0 


(12) 
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This  is  now  in  a  form  suitable  for  application  of  the  WKB  technique.  One  obtains  approx¬ 
imately: 


sin  /  [(1  +  Z )"'  +  (1  +  Z)-2(-KZ'/2  -K2/4)]’'2  dT 
[(1  +  Z)''  +  (1  +  Z)  2  (-KZ’/2  -K2/4)] ,M 


(13) 


Since  both  K  and  Z  are  smaller  than  unity  in  practice  and  as  the  denominator  in  (13)  is 
raised  to  a  small  power,  we  retain  only  the  first  term  therein.  Hence,  (9)  can  be  written 
in  terms  of  (13): 


d2Z 

dT2 


P2 

2(1  +  Z)3/2 


(1  -cos2U).exp.  (-K 


dT 


f 

M  + 


(14) 


where  U  is  the  argument  of  the  sine  in  (13).  To  determine  the  final  velocity  of  the  cone, 
(14)  will  be  integrated  from  0  to  infinity.  The  oscillatory  term  contributes  very  little 
and  can  be  neglected.  Furthermore,  in  accordance  with  the  third  and  fourth  initial  con¬ 
ditions  (10)  we  assume: 


Z  =  q2T2  (15) 

where  q  is  an  as  yet  undetermined  constant.  Changing  the  independent  variable  as 
follows: 


qT  =  tan  6 


we  obtain  for  the  integral  of  (14): 


<Z’) 

CO 


cos6.exp.  (-kq'1  6).d<5 


which  can  be  evaluated  to  give: 

(Z,)«  =  [K/q  +  exp.(-K7r/2q)]  .(1  +  K2/qV 


(16) 


from  which  the  transfer  efficiency  can  be  determined  as  follows: 


(Z')2/P2  =  2e  =  .  [  K/q  +  exp.(-K*r /2q)]*  .  (1  +  KJ/qV* 


(17) 
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The  constant  q  is  now  determined  by  considering  that  as  the  loss  factor  K  vanishes,  the 
efficiency  must  approach  unity.  It  follows  that  q  =  P.(8)'U2. 

It  appears  from  (17)  that  the  transfer  efficiency  is  a  function  of  F/K,  as  indeed 
it  must  according  to  (4) .  Moreover,  comparison  of  the  exact  expression  (4)  with  the 
approximation  (17)  in  Figure  1  shows  that  agreement  is  good  over  most  of  the  range  of  the 
function.  Also  shown  in  Figure  1  are  the  points  obtained  from  the  solution  of  the  system 
(8)  -  (10)  by  means  of  an  analog  computer.  The  agreement  is  good  except  for  low  values 
of  K  where  the  available  integration  time  is  insufficient  to  yield  a  final  constant  value  for 
the  efficiency. 

Further  comparison  between  the  analog  and  the  WKB  solutions  was  attempted 
by  determining  the  value  of  T  for  which  the  current  Y  passes  through  zero.  This  value 
was  determined  graphically  from  the  analog  solution  and  compared  with  the  WKB  solution 
(13)  and  (15).  The  comparison  is  shown  in  Figure  2  where  it  can  be  seen  that  the  agree¬ 
ment  is  adequate  for  the  greater  part  of  the  range  of  the  parameters  and  variables  involved. 

It  follows  from  Equation  (4)  and  Figure  lthat  the  efficiency  can  be  kept  as  close 
to  unity  as  desired  provided  the  condenser  voltage  is  chosen  high  enough.  Physically  this 
may  be  interpreted  to  mean,  that  the  velocity  with  which  the  cone  leaves  the  coil  deter¬ 
mines  the  fraction  of  the  total  energy  which  goes  into  ohmio  loss;  if  this  velocity  is  high, 
the  time  available  for  creation  of  this  loss  is  short  compared  with  the  time  taken  to  accel¬ 
erate  the  cone  to  its  final  speed. 

A  much  more  serious  limitation  on  the  efficiency  attainable  is  imposed  by  the 
fact  that  in  practice  it  is  difficult  to  achieve  a  configuration  of  electromagnetic  acceler¬ 
ator  in  which  the  inductance  is  a  strictly  linear  function  of  the  load  displacement.  In  the 
present  case,  measurement  of  the  variable  inductance  with  the  aid  of  an  impedance  bridge 
showed  a  functional  dependence  closely  following  that  of  a  negative  exponential  (Fig.  3) . 

We  shall  show  in  the  following  to  what  extent  the  previous  analysis  must  be  modified  to 
take  account  of  the  non-linear  behaviour  of  the  circuit  inductance. 

2. 1  Non-Linear  Change  of  Inductance  with  Load  Displacement 

From  the  above  mentioned  inductance  measurements  it  appears  that  the  depend¬ 
ence  of  the  total  circuit  inductance  on  the  load  displacement  can  be  represented  by: 

Li 

L  =  L,  +  —  [  1  -  exp.(-dx)  ] 


or  in  terms  of  the  non-dimensional  parameters  (Q): 


L  =  L,  [l  +  Q-'  (1  -  exp.(-QZ))] 


(18) 


where 


Q  =  L0d/L, 


(19) 
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is  the  non-dimensional  inductance  coefficient.  Putting: 


1  +  Q  '[l -exp.(-QZ)]  =  f(Z) 


we  note  that  (8)  remains  unchanged  except  that  f(Z)  is  substituted  for  1  +  Z.  In  (9),  how¬ 
ever,  we  must  replace  L,  by  the  more  general  term  dL/dx,  whence  we  obtain: 

0  =  Y 2 .  exp .  (-  Q  Z)  (20) 


Proceeding  as  before  by  utilizing  the  WKB  technique,  we  obtain: 

—  =  2-'P2r3,2exp.(-QZ  -K  /  f-'dT)  (21) 

dT 

From  this,  the  maximum  transfer  efficiency  emoi(  can  be  derived  at  once.  We  note  first 
that  as  K  =  0,  the  differential  Equation  (21)  does  not  contain  T  explicitly;  one  may 
therefore  put: 

dZ  .  „  .  d Ei 

dT  p’  dT’  1/2  dT 

and  hence  the  limiting  velocity  pm  is  given  by: 


p2  =  -  p2Q-'  f  (l  +  Q-'  -  Q-'exp.(-QZ))'3,2d(exp.-QZ) 

rp  ^ 

o 


from  which  follows  that: 


P"2p2 

*  m 


=  2c, 


=  2(1 


4 


R. 


Q  +  1 


(22) 


It  will  be  seen  from  (22)  that,  contrary  to  the  linear  case,  the  maximum  efficiency  can  no 
longer  be  made  as  close  to  unity  as  desired.  Figure  4  shows  that  the  efficiency  drops 
away  from  unity  to  zero  as  Q  varies  from  0  to  infinity. 

Of  greater  practical  interest  is  the  general  case  where  K  does  not  vanish.  An 
approximate  solution  may  be  obtained  using  the  same  variable  substitution  as  in  2.0.  We 
derive: 

f  *;2 

e  =  J  cos  6. 


exp.(-Kq''6  -Q  tan2 6) 


(23) 
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where  q  is  again  given  by  P.  (8)'1'2.  This  integral  may  be  evaluated  explicitely  by  approx 
imating  exp.  (-Q. tan 2 6)  by  an  even-term  Fourier  series  of  three  terms,  thus: 


exp.  (-Q.tanJ6)  =  a0  +  a2cos  26  +  a4cos  4  6 


where 


»/  2 


=  4 . 7T  ' 


/  cos  n6.exp. (-Q.tan56)  d6 


It  is  convenient  to  transform  these  integrals  into  those  involving  only  powers  of  the  cosine 

th  i 

»/2  «. 
bn  =  /  cos2n6.exp.(-Qtan26).d6  =  1/2  /  u-1'2  (1  +  u)-n'  exp.(-Qu)  du 

o  o 


The  last  integral  on  the  right  may  be  expressed  in  terms  of  confluent  hypergeometric 
functions  which  in  this  case  reduce'01  to  combinations  of  square  roots,  exponentials  and 
error  functions.  Finally  it  remains  to  evaluate  integrals  involving  products  of  exponen¬ 
tials  and  cosines,  the  result  of  which  is  easily  obtained'". 

Results  of  these  calculations  are  shown  in  Figure  5,  where  it  is  seen  that  the 
efficiency  in  each  curve  reaches  a  saturation  value  which  cannot  be  exceeded  no  matter 
how  large  the  driving  energy.  As  an  additional  check,  the  saturation  efficiencies  were 
plotted  in  Figure  4  and  compared  with  the  simpler  formula  (22) .  It  is  seen  that  although 
the  three-term  Fourier  series  is  a  rather  rough  approximation,  the  agreement  is  fair. 

In  conclusion  of  this  chapter  we  shall  derive  one  additional  formula  which  has 
proved  particularly  useful  in  the  comparison  of  theory  and  experiment.  In  most  of  our 
experiments  the  ratio  P/K  was  quite  large.  The  first  exponential  in  (23)  will  therefore 
reduce  to  small  values  for  small  6.  It  is  therefore  admissible  to  replace  the  ine  by 
unity  and  tan  6  by  62,  obtaining: 


f. 


exp.(-K.q-' 6  -Q6 


2 


The  same  restrictions  imposed  on  the  variables  also  permits  us  to  replace  the  upper  in¬ 
tegration  limit  by  infinity.  The  integral  can  be  evaluated  using  known  formulas'2 '  yield¬ 
ing: 


e  ~  n  (4Q)''  ^exp.^'j^rfc^-''2)]2 


where: 


w  =  P2Q/2K2 
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Further,  if  w  is  small,  the  asymptotic  limit  of  the  error  function  complement  may  be 
used,  giving: 


e  =  —  (1  -  Qw) 


(24) 


In  Figure  G  experimental  efficiencies  were  compared  with  calculated  ones  using  (24).  The 
agreement  is  seen  to  be  quite  good.  The  experimental  techniques  used  in  obtaining  the 
measured  values  will  be  discussed  below. 

In  the  foregoing  it  has  been  assumed  that  the  mass  of  the  cone  plus  the  virtual 
mass  of  the  liquid  contained  are  independent  of  the  speed  and  acceleration,  which  is 
generally  the  case. 


3.0  HYDRODYNAMICS  OF  JET  FORMATION 
3. 1  Stability  of  Rotating  Liquid  in  a  Cone 

If  a  circular  cone  partially  filled  with  liquid  is  rotated  about  its  axis,  the 
meniscus  of  the  liquid  rises  to  form  a  paraboloidal  free  surface.  Contrary  to  what 
occurs  in  a  vessel  with  cylindrical  sides,  this  process  cannot  be  continued  to  ever- 
increasing  rotational  speeds.  When  this  speed  reaches  a  value  at  which  the  paraboloid 
becomes  tangent  to  the  side  of  the  cone,  the  liquid  becomes  unstable  and  spills  over  the 
edge  of  the  cone.  Thus  the  limit  to  which  a  liquid  layer  can  be  created  in  this  manner  is 
imposed  by  the  existence  of  a  critical  rotational  speed .  We  shall  calculate  this  speed 
from  simple  geometrical  considerations. 

Figure  7  shows  the  cross  section  of  a  half-cone  with  apex  angle  80  .  The  height 
of  the  liquid  at  zero  rotational  speed  is  H0  and  the  radius  of  its  meniscus  R0.  When 
rotated  it  rises  to  a  height  of  H,  and  a  radius  of  R, .  Pressure  balance  at  the  parab¬ 
oloidal  meniscus  demands  that  the  meniscus  assumes  a  shape  given  by: 

g(H,  -  r,  -  z)  =  1/2 n2(y2  -  rJ) 

Volume  integration  shows  that  the  volume  of  the  liquid  is: 

kJitH,  R  JirB* 

3  4g 


where  ft  is  the  angular  velocity  and  g  the  acceleration  of  gravity.  This  volume  must  equal 
the  volume  of  the  unrotated  liquid: 


R2*H 

O  0 


3 
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In  addition  we  have  from  geometrical  considerations: 


tan  60  =  R,/H,  -  R0/H0 


From  this  and  the  volume  balance  equation  we  obtain: 


R?  -  4gR*  cot  0o/3fi*  +  4R3g  cot  6J 3fi2  =  0 


This  quartic  equation  in  R,  has  the  form: 


X4  -  MX3  +  N  =  0 


The  meniscus  is  tangent  to  the  cone  provided  that  this  equation  has  two  coinciding  roots, 
the  condition  for  which  is: 


27M4  =  256N 


from  which  is  obtained  the  critical  rotational  speed  fi*: 


n*  =  2'u3cot  9  .  (g/H0)'12 


(25) 


This  relation  is  shown  graphically  in  Figure  8  where  the  critical  speed  is  plotted  against 
the  initial  depth  h  =  H,  -  He.  Also  shown  are  experimental  points  obtained  with  water 
and  mercury.  The  agreement  may  be  considered  fair  in  view  of  the  presence  of  capillary 
effects  and  the  difficulty  in  determining  the  exact  point  of  criticality. 

3.2  Jet  Formation  Under  Artificial  Gravity 

We  shall  here  diverge  somewhat  from  the  classical  theory  of  Monroe  jet  for¬ 
mation,  as  in  our  case  the  liquid  layer  is  not  of  uniform  thickness.  An  observer  moving 
with  the  cone  will  conclude  that  the  meniscus  is  the  proper  equilibrium  shape  between 
the  actions  of  gravitation  and  centrifugal  force.  As  the  cone  accelerates  upward,  the 
observer  will  experience  a  sudden  increase  of  apparent  gravity  and  will  rightly  conclude 
that  the  meniscus  must  assume  a  different  shape.  The  unsteady  motion  of  this  readjust¬ 
ment  will  cause  a  jet  to  develop  near  the  trough  of  the  meniscus  which  will  rapidly  move 
ahead  of  the  edge  of  the  cone.  How  this  motion  is  influenced  by  the  rotational  speed  and 
the  filling  level  of  the  cone  will  be  analyzed  in  the  following. 

Since  for  practical  reasons  the  interest  centers  mainly  on  a  knowledge  of  the  ve¬ 
locity  of  the  head  of  the  jet,  much  valuable  information  can  be  gathered  from  a  linearized 
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analysis.  The  full  problem  involving  the  non-linear  terms  in  the  equations  of  hydrodyna¬ 
mics  as  well  as  considerations  of  a  free  surface  of  ever-changing  shape  requires  extensive 
numerical  study  and  will  therefore  not  be  considered  here.  On  the  other  hand  Penney  and 
Thornhill'31  have  shown  that  a  linearized  treatment  with  initially  fixed  surfaces  may  yield 
valuable  information  relating  to  such  problems  as  the  collapse  of  a  fluid  column  under 
gravity.  The  approach  will  be  adopted  in  the  following.  The  fluid  in  the  cone  is  consid¬ 
ered  ideal,  inviscid  and  incompressible  and,  while  the  paraboloidal  meniscus  is  formed 
by  rotational  effects,  all  further  influence  of  rotations  on  th*  motion  is  neglected.  This 
latter  assumption  may  be  subject  to  revision  in  certain  cas  is,  as  will  be  shown  below. 

The  essence  of  the  method  of  Penney  et.al.  lies  in  the  neglect  of  the  velocity 
convection  terms  in  Euler'  equation,  toget\er  with  considering  the  free  surface  to  be 
acted  upon  by  the  hydrostatic  pressure  of  the  fluid  column  above  the  lowest  point.  The 
fluid  adjacent  to  the  other  bounding  surfaces  are  subject  to  the  usual  conditions  of  van¬ 
ishing  transverse  flow.  If  all  surfaces  involved  are  members  of  the  same  orthogonal  set 
of  co-ordinate  surfaces,  the  problem  is  solved  at  once  by  the  method  of  eigenfunctions 
expansion.  Though  unfortunately  the  paraboloidal  meniscus  and  the  conical  container  do 
not  belong  to  the  same  co-ordinate  system,  a  two-step  approximation  procedure  has  been 
adopted  which  was  found  to  give  useful  information  of  the  motion  for  the  case  where  the 
distortion  of  the  free  surface  is  not  too  great.  However,  to  permit  this  last  step,  it  is 
necessary  to  make  the  additional  assumption  that  the  cone  is  slender,  i.e.  that  its  apex 
angle  is  small. 

3. 2. 1  First  Order  Approximation 

We  adopt  a  spherical  polar  co-ordinate  system  as  depicted  in  Figure  7.  We 
seek  a  solution  of  Laplace's  equation  for  the  velocity  potential  <p  subject  to  the  boundary 
conditions: 


8  <p 
80 


0 


for  0 


8<p 

at 


free  surface 


p/p  =  z.dV/dt 


(26) 


The  first  of  these  conditions  expresses  the  fact  that  no  flow  exists  across  the  surface  of 
the  cone.  The  second,  in  which  p  is  the  pressure,  p  the  density  and  V  the  vertical  veloc¬ 
ity  of  the  cone  derives  from  the  condition  of  Penney  et.al.  that  the  pressure  on  the  free 
surface  is  equal  the  hydrostatic  pressure  due  to  the  height  of  the  fluid  above  the  parabolic 
trough*.  In  the  first  order,  the  free  surface  is  defined  as  r  =  r,.  This  is  in  accordance 
with  the  assumption  that  the  cone  is  slender  and  that  the  constant  radius  vector  describes 
an  essentially  plane  surface.  The  general  solution  of  Laplace's  equation  can  be  written 
thus: 


<p  =  Tj  An(r/R,)".Pn(cos  9) 


*  It  is  assumed  that  dV/dt  is  much  larger  than  g. 
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where  An  are  undetermined  constants  and  Pn  are  Legendre  functions.  Since  0o  is  assumed 
small,  all  the  running  indices  (which  may  be  non-integer)  will  be  large.  It  is  therefore 
permissible  to  substitute  the  asymptotic  Bessel  function  representation  for  the  Legendre 
functions,  thus: 


Pn  (cos  0)  ^  Jo  (n  + 


to  obtain: 


<P  =  £An(r/H,)q"/0°"1/2  J(qn0/0o) 


(27) 


where  qn  are  the  zeros  of  the  Bessel  function  J,(q„)  =  0.  Writing  the  equation  of  the 
parabolic  surface  in  Figure  7  as: 


z  =  yy2  where  y  =  fl2/2g 


we  see  that: 


(z/y)'12  =  (r,  +  z).tan  9  £  (r,  +  z ).B 


where  the  tan  has  been  replaced  by  its  argument  in  accordance  with  the  slender  cone 
assumption.  Solving  for  z  to  second  order: 


z  £  y02r2 


obtaining  for  the  second  condition  (26): 


free  surface 


y  02r2.dV/dt 


Integrating  with  respect  to  t, 


<P 


r 


R, 


=  yV#2r2 


(28) 
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This  is  the  final  form  of  the  second  boundary  condition  (26).  Together  with  (27)  we  there¬ 
fore  have: 


Z  An.Jo(qn0./0o)  =  yV0lrl(6/eo)7 


This  is  Dini's  form  of  a  Fourier-Bessel  expansion  which  can  be  solved  immediately  by 
employing  the  series  representation: 


Z  JotanXM'oV.Jofan)  =  X  J/4 


from  which  follows  that  the  coefficients  of  the  series  are: 


An  =  4Vy^rjq-nVJ0(qn) 


To  determine  the  velocity  components  on  the  free  surface,  we  remember  that: 


vr  =  d<p/dr,  vs  =  d<p/(r.d6) 


whence  on  the  free  surface: 


vr  =  4V  y  0„r, .  F(0/8„) 


v9  =  2V  y  0or}(0/0o) 


(29) 


The  function  F(x)  is  the  sum  of  the  series: 


F<x)  =  Z  Jo^nXi-q;' /J0(qn) 


(30) 


where  the  term  -1/2  has  been  neglected  relative  to  the  much  larger  values  of  qn/0o. 


3.2.2  Approximate  Summation  of  the  Series  F(x) 

With  (29)  the  problem  of  the  first-order  solution  is  formally  solved;  however, 
the  series  (30)  is  cumbersome  and  converges  slowly.  We  shall  attempt  an  approximate 
summation  of  this  series  in  finite  terms.  Since  q„  are  large  numbers,  it  is  admissible 
to  represent  J0(q„)  in  terms  of  its  asymptotic  value: 


J0(q„)  ^  (2/irqn),t2  cos  (qn  -  ir/4) 


(31) 


-  14  - 


and  as  approximately: 


q„  =  *(n  +  1/4) 

it  follows  that  (30)  may  be  written  as: 

F(x)  =  (*/2)  Yj  Jo(q„x).J0(qn) 

On  the  other  hand  we  have: 

J0(xn)  ~  (x)*' (2/n)IIJ  cos  (ffn  -  tr/4)  =  n‘W2/ff 


substituting  in  (31)  for  q„  =  ff  (n  +  1/4)  it  appears  that  on  neglecting  1/4  under  the 
square  root,  there  results: 


J0(irn)  £  2-"2J(qn)  (32) 


If  in  (30)  the  argument  q„x  is  replaced  by  x  *(n  +  1/4)  and  the  addition  theorem  of  the 
Bessel  functions ,4)  to  the  first  two  terms  is  applied,  we  obtain: 

F(x)  =  (ir/2"2)  Y  [jQ(ffn)- Jo(xirn).  J0(l/4xff)  -2U2  J0(qn)J,  (xirn)J,(l/4x)] 


Consider  the  first  term  in  this  expression,  i.e.: 

F,  (x)  =  (jt/2)  Y  J8(xJrn).J0(irn) 


Applying  the  product  representation  of  Bessel  functions 15 ’: 


F,(x)  =  1/2  Y  f  J0 +  xJ-2x.cos  u)1'2 J  .du 
and  the  Schloemilch  series  ,4): 

Y  J0(nx)  =  -1/2  +  1/x 
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we  can  sum  the  infinite  series  before  evaluating  the  integral,  further  we  use  the  definite 
integral 1 7): 

n 

f  (1  -  2p.oos  u  +  p2)"',2du  =  2K(p) 

o 

where  K(p)  is  the  elliptic  integral  of  the  first  kind.  Hence  we  find: 

F,  (x)  =  1/2  [(-1/2  )»  +  (2/tt)  K  (x)  ]  (33) 

Since  the  remaining  term  is  a  first-order  correction,  we  are  free  to  replace  xn  by  qn, 
neglecting  again  x/4,  we  get: 

Fj(X)  =  Z  JofanMlton) 
or,  using  the  asymptotic  value  for  Jo(qn)  again: 

Fj(x)  =  Z  (2A)  J,(xqn).q;'/Je(q„) 

This  is  a  summable  Dini  series  whose  result  is: 


F2(x)  =  -x  !x 


Substituting  F,  (x)  and  F2(x)  in  the  original  expression,  we  get: 

F(x)  =  2',(2  J0(xx/4).  [— jt 2 _1  +  (2/*)K(x)]  +  x.J,(xx/4)  (34) 


As  a  numerical  check  we  calculate  F(0)  =  .40361  and  compare  this  with  the  value  ob¬ 
tained  by  direct  summation  of  the  series,  giving  .38485.  The  error  is  less  than  5%. 

3.2.3  Virtual  Mass  Ratio 

The  final  efficiency  departs  from  the  transfer  efficiency  calculated  above  by  the 
fact  that  along  with  the  fluid,  a  metallic  cone  has  to  be  accelerated,  which  by  itself  does 
not  contribute  to  the  penetrating  action  of  the  jet.  Moreover,  the  mass  m  in  (4)  is  not  to 
be  taken  as  the  actual  rest  mass  of  the  cone-fluid  system  but  must  contain  the  virtual 
mass  of  the  fluid  created  by  the  flow  itself.  We  shall  now  derive  a  first-order  expression 
for  the  ratio  of  virtual  to  actual  mass  of  the  fluid. 
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The  vertical  component  of  the  integrated  pressure  forces  acting  on  the  side  of 
the  cone  is: 


f,  =  27T 


/  rp.  sinJ0odr 


o 


and  from  the  expression  for  the  pressure  in  terms  of  the  velocity  potential  and  the  solu¬ 
tion  (29): 


f.  = 


87ry0  Vp 


dT  E  Or  + 


The  mass  force  is: 


f 


m 


m 


dV 

‘dt 


whence  the  virtual  mass  ratio: 


f,  +  tm 


1  +  12S2Jr1g'1  Yj  q;2<r  +  3/2)-1 


The  infinite  series  may  be  summed  in  terms  of  the  if>  function  (logarithmic  derivative  of 
the  gamma  function)  to  give: 


f  +  f 

— - — -  =  1  +  «2r,g-'  [l  -  48(9tt0o)-  1  (*(5/4  +  3»./2ir)  -  ft (5/4)1  (35) 

Substituting  (25)  for  the  angular  velocity,  one  sees  that  in  the  case  of  a  12 “cone,  the 
virtual  mass  is  increased  by  more  than  60%. 

3.3  Second  Order  Approximation 

The  first  order  approximation  was  derived  on  the  assumption  that  only  the 
pressure  force  was  influenced  by  the  curvature  of  the  free  surface,  but  that  the  resulting 
flow  proceeded  as  if  the  interface  were  plane.  We  shall  now  attempt  to  improve  upon 
this  approximation  by  partially  taking  account  of  the  fact  that  the  pressure  condition  (26) 
is  to  be  applied  to  the  actual  shape  of  the  meniscus  rather  than  to  its  rest  condition.  We 
write  the  equation  for  the  parabolic  meniscus  in  polar  co-ordinates  as  follows: 

r/H0  =  1  -  a  +  0  (0  2/e\) 


(36) 
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where  a  and  (3  are  constants  depending  on  the  speed  of  rotation.  In  the  following  we  shall 
utilize  a  non-dimensional  parameter  S  which  fully  describes  the  motion,  together  with 
the  cone  apex  angle  0O: 


S  =  02flo  H0/(2g)  =  Fr/(8RoJ) 


(37) 


where  Fr  is  the  Froude  number  and  Ro  the  Rossby  number  referred  to  the  characteristic 
length  0OHO.  From  the  geometrical  considerations  of  paragraph  3. 1  it  follows  that: 


a  =  seo/2,  p  =  se0 


(38) 


The  procedure  to  be  followed  subsequently  involves  the  substitution  of  (36)  in  (27),  where 
we  have  adopted  Hcas  reference  length  in  place  of  H,,  expanding  the  resulting  trinomial 
to  first  order  according  to  the  binomial  theorem  and  replacing  A„  by  the  sum  of  the  two 
coefficients  representing  the  coefficients  of  the  Dini  series  for  x2  and  x4  multiplied  by 
unknown  constants  a  and  b  respectively.  Since  the  whole  series  must  equal  the  pressure 
on  the  meniscus  and  therefore  be  proportional  to  x2,  we  must  adjust  the  constants  a  and 
b  in  such  a  way  that  any  terms  in  x4  must  vanish  and  that  the  coefficient  of  x2  correspond 
to  that  on  the  right  hand  side.  Neglecting  -1/2  again  relative  to  qn/0o,  expanding  (36) 
raised  to  the  power  qn/0o  to  two  terms  and  substituting  in  (27),  equating  this  to  the  pres¬ 
sure  (28),  we  obtain: 


I  (a"  +  b")  [l-(qm/0n)(a  -  0X2)J  j0(qnx)  =  VyOl H„xJ 


where: 


a„  =  4aq-n2/J0(qn),  b„  =  8bq;2(l-8/q2)/Jo(qn) 


are  respectively  the  coefficients  of  the  Dini  series  representing  ax2  and  bx4.  Summing 
the  resulting  series,  there  results: 


ax2  —  (or  - 0x2). F(x)  +  bx4  =  Vy6JH2x2 
“0  0  0 


(39) 


where  F(x)  is  the  function  given  in  (34).  Note  that  the  term  resulting  from  multiplication 
of  the  series  for  bx4  with  q„  has  been  omitted,  in  accordance  with  order -of -magnitude 
consideration.  F(x)  must  now  be  expanded  in  terms  of  x2  and  a  and  b  must  be  chosen 
such  that  the  x4  term  disappears  and  the  x2  term  balances  against  the  right  hand  side  of 
(39).  From  (34)  we  find: 


F(x)  =  2'1 2 -,,2(2  -  ir)  +  (l/2)  x2(1  +  ir/4) 
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whence  the  coefficients  a  and  b: 

V  y^HJ 

'  o  o 

a  -  - — - — - 

1  +  2/3(0o)-'2-’/2(2  -  ?r)  -  2a0;,(l  4  tt/4) 


2/3  Vy0oH2(l  +  4 2 1/2  +  4.2U2  (?r  -  2)) 

b  =  - - 

1  +  2/3fl;l2‘,,2(2  -  IT)  -  2o0;’(l  4  TT/4) 


Now,  the  quantity  of  greatest  interest  is  the  velocity  of  the  head  of  the  jet,  i.e. 

vf0  =  vr  (®  ~  0)  ~  &$/ 3  r  (0  =  0).  Returning  to  the  series  representation  of  the  velocity 
potential: 


v 


'O 


y^a-ai'1-8"-372. 


+ 


8b 

0oQnJO(qn) 


(39a) 


Note  that  in  the  second  term  the  1  -  a  term  has  been  omitted  to  keep  orders  of  magnitude 
correct.  The  first  term  in  this  may  be  written,  using  the  approximate  value  for  q  and 
the  asymptotic  expression  for  J0(q„): 


I  (ff/2).exp.[(7rnr'  4  (i/4)tt 0 * ’ -3/2)log(l  -  o)j  . Jo(7rn). 2W2 
We  expand  log  (1  -  or)  =  -  a  and  use  the  Schloemilch  series 181  to  the  second  order: 

I  J0(nir).exp:(nu)  =  -1/2  4-  (UJ  +  w2)’'12  +  u/12 


The  second  term  in  (39a)  was  summed  term  by  term  and  gave  the  value  .  77021.  b.  Taking 
all  terms  together  and  expanding  a  and  b  to  first  order  in  a  and  /3  we  get  finally: 

vfo/V  -  2. 2,/2  (-1  4  tt/2)S  (1  -  SU>  (40) 


where: 


U 


24  (-1  +  7T/2) 


(»/ 8) -(3/4)0o -l-(7T/4)  -2 ' ' 2  (IT  -  2)  4 


.77021  (1  4  tt/4) 
2,,J  (-1  4  tt/2) 


This  is  the  second-order  expression  for  the  velocity  amplifeation  ratio,  i.e.  the  factor 
by  which  the  velocity  of  the  cone  is  multiplied  to  obtain  the  velocity  of  the  head  of  the  jet. 
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It  must  be  remembered  that  the  whole  treatment  was  carried  out  in  a  co-ordinate  system 
travelling  with  the  cone.  To  convert  to  a  system  at  rest  in  the  laboratory  frame,  one 
must  add  unity  to  (40) .  The  largest  possible  value  of  this  ratio  is  obtained  for  the  case  of 
critical  rotational  velocity  (25).  In  this  case  S  =  Sm0,  =  . 31498/0O.  In  our  experiments 
#o  -  tt/15.  This  yields  in  the  laboratory  frame: 


(v,0/V)mo,  -  6.12958 


It  is  interesting  to  note  that  in  first  order  of  Sma>  vfo/V  becomes  inversely  proportional 
to  60  this  is  in  agreement  with  the  findings  of  the  classical  theory  of  shaped  charges 

3.4  Shape  of  Jet  Profile  for  Small  Times 

From  (29)  it  is  possible  to  obtain  the  shape  of  the  jet  by  simply  integrating 
these  equations  with  respect  to  time.  It  must  be  remembered  that  in  this  treatment  the 
convection  terms  in  Euler's  equation  were  neglected;  all  solutions  apply  therefore  only  to 
small  times.  Hence,  it  is  not  necessary  to  proceed  to  the  second  order  in  calculating  the 
jet  profile,  as  the  first  order  will  yield  enough  information  for  the  case  of  low  distortion 
of  the  free  surface. 

Figure  9  shows  a  number  of  profiles  calculated  for  progressively  larger  values 
of  time.  A  distinctive  feature  of  these  curves  is  the  rapid  negative  motion  of  the  fluid 
near  the  cone  wall.  This  is  due  to  the  logarithmic  infinity  of  K(x)  in  (34)  as  x  approaches 
one.  In  the  experiments  some  breakup  of  the  fluid  is  observed  near  this  region,  probably 
due  to  cavitation  caused  by  the  large  negative  pressures  generated.  In  Figure  9  are  also 
shown  some  experimental  points  taken  from  photographs  of  the  developing  jet.  The  agree¬ 
ment  appears  to  be  satisfactory  in  view  of  the  assumptions  made  and  the  difficulty  in 
maintaining  good  experimental  conditions.  Poorest  agreement  is  found  at  earliest  times, 
probably  due  to  the  fact  that  the  initial  free  surface  was  not  plane,  as  assumed  in  the 
first-order  theory,  but  distorted  by  capillary  action  near  the  cone  wall. 

3.5  Stream  Lines 

The  stream  lines  of  the  first-order  solution  can  be  obtained  from  the  definition 
of  the  stream  function  i p: 


\  =  r-*.0-'dil>/d0 

v,  =  -(rO)  ' dip/dr 


hence  from  (29): 


(qn/0e)  +  i/2 


<P  =  E  A„  H,«(l-«./2qn).(r/H,) 


•J.(q„V0„) 


(41) 
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where  An  are  the  coefficients  of  the  Dini  series  given  before.  The  stream  lines  arc  shown 
plotted  in  Figure  10  for  the  case  of  a  cone  of  24°  included  apex  angle  ( 60  =  12°).  The 
large  negative  velocities  near  the  cone  well  are  again  apparent  in  this  figure,  as  well  as 
the  cumulation  of  stream  lines  near  the  axis  of  the  cone  where  the  jet  is  formed.  It  must 
be  noted  that  the  stream  function  (41)  applies  to  a  frame  of  reference  fixed  in  the  moving 
cone;  for  a  frame  stationary  in  the  laboratory,  the  term: 


ipQ  =  -Vr20  2/2 


must  be  added  to  (41). 


4.0  EXPERIMENTS 

4. 1  General  Considerations 

The  experimental  rig  consists  essentially  of  a  condenser  bank  discharging  into 
a  conical  single-turn  coil  carrying  the  rotating  cone.  The  similarity  with  a  plasma  pinch 
(theta  pinch)  experiment  is  obvious;  however  the  design  criteria  are  somewhat  different. 
As  in  the  plasma  pinch,  the  time  constant  of  the  discharge  must  be  kept  small  to  ensure 
that  the  magnetic  field  of  the  coil  does  not  penetrate  through  the  wall  of  the  cone.  This  is 
usually  achieved  by  keeping  the  circuit  inductance  Leas  small  as  possible,  a  requirement 
which  also  guarantees  that  the  inductance  coefficient  Q  (19)  remains  small.  On  the  other 
hand  it  is  seen  from  (4)  that  for  greater  transfer  efficiency  the  variable  inductance  L, 
must  be  made  large.  Since  most  of  the  circuit  irductance  is  located  in  the  condenser 
bank  and  the  external  conductors,  a  transformer  designed  to  match  Li  into  the  primary 
circuit  could  be  used  to  optimize  the  efficiency.  Unfortunately,  the  high  frequency  of  the 
discharge  transient  precludes  the  use  of  iron  transformer  cores;  the  use  of  ferrite  as  core 
material  is  the  next  best  choice,  though  even  in  this  case  the  transformer  develops  a 
considerable  leakage  inductance.  The  compromise  finally  adopted  comprised  a  ferrite- 
cored  3tep-down  transformer  of  ratio  40:1. 

4.2  Transformer  and  Coil 

The  schematic  of  the  transformer  is  shown  in  Figure  11.  The  ferrite  core  is 
eveloped  in  a  thin  layer  of  mylar  sheet  insulation  which  carries  the  40-turn  primary 
wound  with  square-section  magnet  wire.  A  further  layer  of  insulation  separates  this 
from  the  secondary  which  is  a  1/8  inch  thick  split  copper  tube,  which  in  turn  is  encased 
in  a  1/2  inch  thick  cylinder  of  steel  to  prevent  the  coil  from  spreading  under  the  action  of 
the  magnetic  pressure.  The  conical  pinch  coil  is  silver-soldered  to  the  gap  in  the  second¬ 
ary  and  is  held  firmly  to  the  transformer  by  means  of  a  1/2  inch  retainer  plate  bolted  to 
two  jaws  welded  to  the  steel  cylinder.  The  coil  itself  is  a  split  brass  block  of  square 
section,  carrying  a  conical  hole  into  which  fits  the  rotating  cone.  The  coil  is  fitted  with 
an  air  bearing  supplied  through  a  nozzle,  while  a  second  nozzle  supplies  air  tangentially 
to  the  cone  to  provide  the  torque  necessary  to  set  the  cone  into  rotating  motion.  The 
transformer  is  held  between  two  plexiglass  rings  and  three  steel  bolts.  Figure  12  shows 
a  photograph  of  the  transformer-coil  assembly. 
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The  variable  inductance  L|  was  measured  by  mounting  the  cone,  made  of  alu¬ 
minum  on  a  traversing  rig,  the  coil  being  clamped  firmly  to  the  base  of  the  rig  and  the 
displacement  of  the  cone  being  measured  with  adial  gauge.  The  inductai  oe  was  measured 
on  a  Hewlett-Packard  R-X  meter  (Schering  Bridge)  using  a  suitable  frequency  of  the  order 
of  megahertz.  The  readings  had  to  be  modified  to  take  account  of  the  lower  working  fre¬ 
quencies  in  the  discharge  circuit,  owing  to  the  greater  penetration  of  the  magnetic  field. 
(Fig.  3). 

4.3  Linear  Velocity  Measurements 

The  transformer-and-coil  assembly  was  installed  in  the  discharge  circuit  and 
the  condenser  bank  charged  to  a  voltage  varying  between  1200  and  7100  V.  The  air  cush¬ 
ion  on  which  the  cone  floats  is,  as  a  rule  si  cient  insu’ation  to  prevent  the  cone  from 
shorting  out  the  coil.  For  greater  safety  high  levels  of  condenser  voltage  the  cone 
external  surface  was  anodized. 

The  "muzzle  velocity"  with  which  the  cone  leaves  the  coil  on  discharge  was 
measured  optically  by  projecting  a  collimated  beam  of  light  horizontally  over  the  cone, 
through  a  45 "glass  splitter  plate  and  back  again  as  a  parallel  beam  about  2  inches  above 
the  first.  The  light  from  this  beam,  together  with  a  portion  of  the  first  beam  separated 
from  it  by  a  half-silvered  mirror,  enters  a  photodiode  the  output  of  which  actuates  an 
electronic  interval  counter.  The  discharge  current  is  simultaneously  monitored  with  the 
aid  of  a  Rogowski  oil  around  one  of  the  conductors,  and  an  oscilloscope.  The  damped 
oscillatory  signal  obtained  is  analyzed  to  obtain  the  frequency  and  hence  the  circuit 
inductance  L0,  while  the  damping  yields  the  total  loss  resistance  R.  As  the  transfer 
efficiency  was  found  to  be  quite  small,  the  approximate  expression  (24)  was  used  for 
comparison  with  experimental  values.  Results  are  shown  in  Figure  6.  The  apparent 
break  in  the  experimental  curve  around  4  kV  is  caused  by  the  fact  that  for  operation  with 
higher  voltages,  the  pressure  supplying  the  air  bearing  had  to  be  increased  to  avoid  flash- 
over  between  coil  and  cone,  the  resulting  increase  in  the  air  gap  giving  a  small  change  in 
the  readings. 

4.4  Photographic  Experiments 

To  compare  the  hydrodynamic  calculations  with  experiment,  it  was  necessary 
to  take  high-speed  photographs  of  the  development  of  the  jet.  The  setup  used  is  shown 
in  the  block  diagram  Figure  13.  The  output  from  the  Rogowski  coil  is  fed  into  an  ampli¬ 
fier  whose  output  actuates  the  start  circuit  of  the  interval  counter.  A  second  output 
signal  is  retarded  by  a  delay  generator  which  supplies  a  stop  signal  to  the  counter.  A 
signal  is  taken  from  one  of  the  decade  circuits  of  the  counter,  amplified  and  fed  into  the 
flash  unit  which  supplies  the  high  voltage  needed  to  trigger  a  xenon  flash  tube'95.  Each 
time  the  decade  cycles  through  its  period,  it  emits  a  pulse  which  normally  goes  to  switch 
on  a  digit  light  on  the  counter.  This  pulse  is  thus  made  to  trigger  the  flash  circuit.  Up 
to  2500  Hz  of  repetition  frequency  is  attainable  in  this  way.  A  Tektronix  time  mark 
generator  is  used  to  calibrate  the  delay  circuit. 

The  burst  of  flashes  (duration  1/2  microsecond)  illuminates  the  object  while 
the  film  of  a  NRC-Dudgeon  drum  camera20'1  is  being  exposed.  This  camera  has  a  hollow 
drum  carrying  the  film  on  the  inside  surface  and  mounted  on  air  bearings  and  driven  by 
an  air  turbine.  Simultaneously,  the  rotational  speed  of  the  cone  is  measured  by  means 
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of  a  marker  painted  on  the  edge  of  the  cone  whose  motion  is  detected  by  a  light  and  photo¬ 
transistor.  The  period  between  successive  passages  of  the  marker  past  the  photo  pickup 
is  measured  on  another  interval  counter.  The  rotational  speed  of  the  drum  camera  is 
monitored  on  an  oscilloscope. 

For  better  photographic  contrast  the  cone  was  partially  filled  with  milk,  in¬ 
stead  of  water.  Typical  results  are  shown  in  Figure  14,  representing  shots  at  different 
condenser  voltages  and  cone  rotational  speeds.  The  voltage  seems  to  have  only  minor 
effect  on  the  shape  of  the  jet,  except  at  high  voltage  when  the  head  of  the  jet  seems  to 
develop  a  plume  of  spray.  This  is  seen  more  clearly  on  Figure  15  which  shows  a  single 
flash  picture  taken  with  a  still  camera.  This  plume  appears  to  be  caused  by  the  drag  in 
the  stagnation  point  boundary  layer  of  the  surrounding  air.  The  resulting  loss  of  fluid 
and  hence  of  apparent  forward  speed  of  the  jet  head  causes  some  difficulty  in  interpreting 
the  results,  as  will  be  seen  in  the  following.  In  addition,  the  high  voltage  shots  also  show 
the  development  of  a  thin  precursor  jet.  It  is  thought  that  this  jet  is  produced  by  the 
elastic  contraction  of  the  cone  under  the  magnetic  driving  field.  In  the  appendix  is  given 
a  mathematical  treatment,  supported  by  experiments,  of  the  magnetic  field  inside  a  con¬ 
ical  coil.  It  is  seen  that  the  field  rises  steeply  from  the  upper  edge  of  the  cone  to  the 
apex,  a  course  which  might  account  for  the  small  thickness  of  the  precursor  jet. 

It  has  been  assumed  in  the  mathematical  analysis  that  the  effect  of  rotation  on 
the  motion  of  the  fluid  (as  opposed  to  the  development  of  the  meniscus)  can  be  neglected. 
This  is  certainly  true  for  the  early  phases  of  the  motion  and  at  high  driving  voltages.  At 
low  voltage  and  at  later  times,  however,  Figure  16  shows  definite  indication  of  rotational 
effects.  These  are  magnified  by  the  fact  that,  as  can  be  seen  from  the  stream  lines 
Figure  10,  the  fluid  starting  near  the  wall  of  the  cone  is  forced  toward  the  centre  to  form 
the  jet.  Due  to  conservation  of  rotary  momentum,  the  angular  speed  of  the  fluid  compos¬ 
ing  the  jet  surface  may  attain  considerable  magnitudes.  This  is  shown  in  Figure  16  where 
the  jet  edge  is  seen  to  be  subject  to  typical  centrifugal  instability,  resulting  in  the  ejec¬ 
tion  of  droplets 2n. 

In  Figure  17  is  seen  a  series  of  exposures  depicting  the  penetration  of  a  paper 
target  by  a  jet  produced  by  a  6  kV  discharge.  By  connecting  the  jet  head  position  by  a 
straight  line,  it  can  be  seen  that  penetration  occurs  instantaneously.  The  precursor  jet 
is  eliminated  in  the  process  and  the  plume  is  somewhat  reduced. 

4.5  Velocity  Amplification:  Comparison 

From  a  technical  poini  of  view,  the  most  important  quantity  is  the  factor  by 
which  the  head-of-the-jet  velocity  is  amplified  relative  to  that  of  the  cone,  i.e.  the 
velocity  amplification.  An  analytical  expression  of  this  quantity  is  given  in  (40).  Com¬ 
parison  with  experimental  findings  was  made  simply  by  measuring  the  respective  slopes 
of  the  cone  and  jet  exposures,  Figure  14.  Absolute  values  of  velocity  are  easily  obtained 
bearing  in  mind  the  fact  that  the  flashing  rate  and  hence  the  time  between  exposures  is 
tightly  synchronized  with  the  cycling  rate  of  the  counter  decade.  The  scale  of  the  pictures 
■  s  obtained  by  exposing  a  measuring  sccle  in  place  of  the  jet.  Values  of  velocity  thus 
obtained  are  found  to  agree  closely  with  those  measured  photoelectrically  in  4.3,  taking 
into  account  the  total  mass  of  cone  plus  fluid,  including  the  virtual  mass.  Comparison 
between  theoretical  and  experimentally  determined  velocity  amplification  vs  rotational 
speed  of  the  cone  is  made  in  Figure  18.  It  will  be  seen  that  the  experimental  points 
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agree  with  the  theoretical  curve  for  low  values  of  angular  velocity  but  begin  to  diverge  for 
higher  values,  and  that  the  agreement  ic  maintained  for  higher  angular  speeds  for  results 
obtained  with  lower  driving  voltages.  This  is  explained  in  terms  of  lateral  loss  of  fluid 
of  the  jet  head  as  a  result  of  the  dispersal  of  fluid  in  the  plume.  Note  that  the  curves  are 
displaced  vertically  so  that  they  pass  through  a  point  somewhat  larger  than  unity  as  re¬ 
quired  by  theory.  This  is  due  to  the  fact  that  even  at  zero  rotational  speed  the  meniscus 
is  somewhat  concave,  owing  to  the  presence  of  interfacial  tension  between  the  fluid  and 
the  cone  material. 


5.0  CONCLUSIONS 

It  is  shown  that  the  mechanics  of  the  explosively  driven  Monroe  jet  (shaped 
c!  rge  effect)  can  be  simulated  by  means  of  a  fluid  in  a  rotating  and  vertically  accelerated 
cone.  The  theory  of  the  electromagnetic  propulsion  of  the  cone  by  means  of  a  condenser 
discharge  can  be  developed  in  terms  of  the  WKB  approximation.  The  hydrodynamic 
theory,  assuming  an  ideal,  inviscid  and  incompressible  fluid  is  carried  out  in  a  two-step 
procedure,  where  the  velocity  amplification  results  from  a  second  order  approximation. 
Good  agreement  is  found  with  theory  for  such  conditions  where  the  velocity  of  the  jet  head 
is  comparatively  low;  at  higher  velocities  the  jet  tends  to  be  dispersed  by  frictional  inter¬ 
action  with  the  surrounding  air,  making  the  apparent  jet  velocity  somewhat  lower  than 
predicted. 


From  the  point  of  view  of  technological  application,  the  coupling  between  con¬ 
denser  energy  and  cone  that  could  be  attained  in  the  present  experiments  was  rather 
poor,  resulting  in  a  maximum  apparent  jet  velocity  of  230  ft/sec  (extrapolating  for  plume 
loss,  this  figure  may  be  increased  by  a  factorof  3).  It  is  proposed  to  redesign  the  electro¬ 
magnetic  propulsion  device  in  order  to  improve  both  the  "circuit  Mach  number "  M  and 
the  inductance  coefficient  Q.  This  should  result  in  a  much  greater  jet  velocity  attainable. 
An  intrinsic  drawback  of  this  device  is  the  necessity  to  operate  it  vertically. 
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FIG.  2 :  POSITION  OF  ZEROS  OF  DISCHARGE  CURRENT.  COMPARISON  OF 
ANALOG  AND  WK8  SOLUTIONS 
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F\G  <?  MAXIMUM  TRANSFER  EFFICIENCY  *a  Q  COMPARISON  OF  DIRECT  SOLUTION  AND 

THREE-  ERM  FOURIER  SERIES 
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FIG. 5  TRANSFER  EFFICIENCY  vs  P/K  FOR  DIFFERENT  VALUES  OF  Q 


6.6:  TRANSFER  EFFICIENCY  vs  CONDENSER  VOLTAGE.  COMP/ 
WITH  PHOTOELECTRIC  SPEED  MEASUREMENTS 


G.8  CRITICAL  ROTATIONAL  PERIOD  OF  POOL  OF  LIQUID  IN  CONE  vs  INITIAL  FILLING  DEPT 

COMPARISON  WITH  EXPERIMENTS 


FIG. 9'-  SHAPE  OF  JET  FOR  EARLY  TIMES.  COMPARISON  WITH 

PHOTOGRAPHS 


FIG  12:  PHOTOGRAPH  OF  TRANSFORMER 
AND  COIL 


FIG.  13  BLOCK  DIAGRAM  OF  PHOTOGRAPHIC  EXPERIMENT 


(a)  V0  *  2 k V,  ROT  SPEED  *  7  35  rps 
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(b)  V0  =  3k  V,  ROT  SPEED  =  8  69  rps 


(c)  V0  *  4 k V,  ROT  SPEED  *  7  46rps 


(d)  V0  *  4k V,  ROT  SPEED  *  5  26rps 


(f )  V0*  6k  V,  ROT  SPEED  =8  20  rps 


FIG. 14  PHOTOGRAPHS  OF  DEVELOPING  JET  IN  CONE 
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FIG.  15  STILL  PHOTOGRAPH  OF  JET  WITH  PLUME 
V0  *  6  k  V  ROT  SPEED,  6.70  r  p  s 
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FIG.  16  DEVELOPMENT  OF  JET  AT  LOW  VOLTAGE  SHOWING  ROTATIONAL  INSTABILITY 
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AXIAL  DISTANCE 
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FIG  19  DISTRIBUTION  OF  AXIAL  MAGNETIC  INDUCTION  IN  CONICAL 
COIL.  COMPARISON  BETWEEN  THEORY  AND  MEASUREMENTS 
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APPENDIX 


To  calculate  the  magnetic  induction  on  the  axis  of  a  conical  coil  due  to  an  applied 
potential  constant  along  its  length,  assume  that  the  current  is  confined  to  a  thin  layer  ad¬ 
jacent  to  the  conical  hole.  This  is  Justified  for  fairly  high  frequencies  where  the  skin 
effect  naturally  restricts  the  depth  of  the  current-carrying  layer.  Calling  the  distance 
of  a  point  on  the  axis  from  the  apex  h,  the  distance  of  the  near  edge  h|  and  that  of  the  far 
end  h2,  we  can  calculate  the  field  induction  along  the  cone  axis  B  by  integrating  over  all 
circular  current  loops  of  thickness  dhc  carrying  a  current  i(h0),  thus: 


B  -  (p/2)  /RJi(h0)(R2  +  (h0-h)V3,2dh0  (42) 

n\ 

where  he  is  the  axial  position  of  the  current  loop,  R  is  its  radius  and  p  is  the  permeability 
of  free  space.  If  the  apex  angle  of  the  cone  is  0O,  we  have: 


tan  eo  =  R/h0 


which  allows  elimination  of  R  from  (42) .  To  estimate  the  functional  dependence  of  i(h0) 
on  h„,  we  assume  that  the  frequency  of  the  applied  potential  is  high  enough  to  make  the 
current  distribution  resistance  dependent.  Hence,  if  the  applied  voltage  is  independent 
of  hD,  we  have  from  Ohm's  law: 


i(h0)  =  k/he 


where  k  is  a  constant.  Substituting  this  in  (42)  we  obtain: 


B 


2hho  +  hV3'2  tanJ0o  dh„ 


which  is  explicttely: 


B  =  (p/2h)k  £(h,  -  h)(h2sec20o  -  2hh,  +  h2)'1'2  -  (h2  -  h)(h2secJ0o  -  2hh2  +  hV'2J(43) 

The  course  of  the  function  B(h),  (43)  is  shown  in  Figure  19.  Comparison  is  made  with 
measurements  using  a  small  magnetic  probe  consisting  of  a  coil  1mm  diameter  and  3  mm 
long  mounted  on  a  traversing  rig.  The  coil  was  energized  by  a  3  Mhz  signal  generator 
matched  into  the  load  by  a  ferrite  cored  step-down  transformer,  The  frequency  of  3  Mhz 
was  chosen  to  obtain  adequate  sensitivity  of  the  probe,  but  variation  of  the  frequency  did 
not  seem  to  cause  material  changes  in  the  results.  As  k  is  an  arbitrary  constant,  the 
theoretical  and  experimental  curves  were  normalized  at  their  peak. 
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As  will  be  seen,  the  agreement  is  good  over  most  of  the  length  of  the  coil.  A 
salient  feature  of  the  result  is  the  rapid  rise  of  the  field  strength  toward  the  apex  of  the 
cone.  When  it  is  remembered  that  the  magnetic  pressure  is  proportional  to  the  square 
of  the  field  strength,  the  dependence  of  pressure  on  axial  position  will  be  a  very  steep 
function.  Calculations  show  that  this  is  one  of  the  conditions  under  which  a  thin  Jet  is 
produced  by  magnetic  pressure  and  elastic  contraction  of  the  cone. 
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